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Abstract: Current theories of nanochannel flow impose no upper bound on flow rates, and predict friction 
through nanochannels can be vanishingly small. We reassess neglecting channel entry effects in extremely 
long channels and find violations at the nanoscale. Even in frictionless nanochannels, end effects provide a 
finite amount of friction. Hence, the speed at which nanochannels transport liquids is limited. Flow-rate and 
slip length measurements are reevaluated using calculations which include end-effect friction. End effects 
are critical for the design of new technological devices and to understand biological transport. 
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Nano s cale channels, such as aquaporin and carbon nanotubes, e xhibit surprisingly large flow rat es ([Preston et al 



19921 : 



1 r 



Hummer et al ,20011 : iKalra et al 



2003 



Falk et al 



20101 : 



Thomas and McGaughey 



20081). Filtration 



membranes fabricated from carbon nanotubes have been reported with flow rates 100 - 100,000 times 



greater than h ave been measured for any other material (Majumder et al 



Du et al 



2005 



20111 : 



Holt et al 



20061 : 



201 1|) . Arbitrarily high measured flow rates are fit to theoretica 



cient 



20061 : 



Y small values for the friction parameter along the channel length (IMajumder et al 



Du et al 



calculations by choosing suffi 



2005 



Holt et al 



201 1|) . These calculations, however, ignore end effects, viscous losses within the liquid near 



the channel's entrance and exit. For nanotubes whose typical lengths are lOOO's times their radii, classical 
fluid dynamics predicts that the great length of the channel provides the primary flow resistance and thus, 
end effects are negligible. The narrow and tortuous channel of aquaporin leads to similar expectations of 
negligible end effects. However, contrary to usual fluid mechanics practice and expectation, our results 
show that long low-friction channels and short channels of any channel friction can have significant end 
effects. Thus, all treatments of nanochannel flow must be aware of end effects. Friction due to end ef- 
fects is essential for calculating physically realistic values for flow rates in nanochannels such as carbon 
nanotubes and aquaporin. Our modified Hagen-Poiseuille equation (see Methods) gives predictions com 



parab 



20091: 



e to molecular dynamics res ults (|Kalra et al 



2003 



Suk and Aluru 



Goldsmith and Martens 



20101 : 



Thomas and McGaughey 



2009|) . In addition, our results provide a theoretical limit to flow rates in any 
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Figure 1. Energy dissipation (high (red), low (purple)), represented by a perceptually balanced color scale (see Methods), 
as numerically computed for a cylindrical channel between two large reservoirs (2-D cross-section, only a small portion 
of the reservoirs is shown). For a traditional high-friction "no-slip" channel (top panel), energy dissipation occurs pre- 
dominantly along the channel length. For a zero-friction channel (bottom panel), energy dissipation occurs solely near the 
entrance and exit of the channel, irrespective of channel length, in regions of high streamline curvature. 

nanochannel, and serve as an upper bound for experimental measurements. 

For short aquaporin or slippery carbon nanotube channels, total resistance within the channel can be 
surprisingly small. Previous efforts to compute flow rates through these channels have focused exclusively 
on the friction within the channel. But, regions of high viscous shear persist near the channel ends, where 
streamline s sharply curve from the large re s ervoirs into (and out of) the small lumen of the channel, see 



Figure 1 (iSampson 



1891 



Weissberg 



1962 



Suk and Aluru 



2010|) . The size of these localized regions is 



small, on the order of the channel radius, yet including them in theoretical calculations shows their contri- 
bution to flow resistance can be significant. In the limit of zero channel friction (/ = 0), flow resistance 
arises solely through shear in these regions near the channel ends, see Figure 1, bottom panel, thus creating 
a fundamental flow-rate speed limit, in contrast to prior theory. 

Figure 2 shows experimental measurements (symbols) of flow rate in channels of various lengths, L. 
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Figure 2. A flow rate speed limit in nanochannels. Predictions (solid lines) of flow rate Q with channel friction / relative 
to the frictio nless-chann el flow rate Qo, versus cha nnel length L, are corn pared with measurements of carb on nanotubes 
(blue circles (iHoltet all 12006) ; red triangles, filled (Maium der et all 120051) and open dMaiumder et al 20111); pu rple stars 
( Du et al 201 ll)) and a quaporin (green square (Preston et al ll992h ). Several reported flow rates dMaiumder e t al 2005J; 
1201 ll; iHolt et a]||2006l) are larger than allo wed for frictionless channels, f = 0. The dashed lin es show predictions, for 
each value of/, when ignoring end effects dMaiumder et aEOOSLlHolt et alll200q ; lDu et all201 ll) . These lines erroneously 
suggest that any Q, no matter how large, is theoretically possible. 



Predictions that ignore end effects (dashed lines) i ncrease without bound , so a valiie of the friction facto r 



/ can be found to fit any empirical measurement (JMajumder et al 



20051: 



Holt et al 



2006 



Du et al 



2011 ) 



Predictions that include end effects (solid lines) asymptote to a maximum flow rate, Qq. The nine mea- 
surements that fall below this line (/ = 0) comply with known physics of fluid flow, while the six above 
this line are unphysically fast. For measurements below the limit, the theory predicts new lower values 
of /, and thus higher values of the slip length (see Methods). For example, new predictions for the two 
overlapping blue circles, see Figure 2, give values for the friction factor (slip length) that are half (twice) as 
large as previously found. For flow measurements closer to the limit Qq, the discrepancy between old and 
new predictions diverges. 

In summary, as friction within channels decreases, end effects make an increasingly large contribution 
to total flow resistance, see Figure 3, even in extremely long channels where end effects have been rotely 
neglected. For short channels, such as aquaporin, end effects remain significant, asymptoting to a finite 
value, regardless of the amount of channel friction, see Figure 3, green line. In previous calculations, which 
omitted end effects, flow rates were unbounded as / — t- 0. Since viscous losses at the channel entry and 
exit are independent of channel friction, including them produces an upper limit to the flow rate (Figure 2). 
Prior measurements that show flow rates above this limit may point to difficulties in accurately determining 
channel radii, net flow rates, or the number of channels spanning the membrane. Errors may also arise due 
to compromise in membrane integrity. Appreciation of end losses is necessary for understanding the high, 
but not unphysically high, flow rates found in carbon nanotubes, so as to ultimately realize the far-reaching 
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Figure 3. Low friction channels are dominated by end effects. As channel friction / decreases, flow resistance R^- due 
to viscous friction at the entrance and exit becomes a larger fraction of the total flow resistance Rt- The colored lines use 
the channel length, channel radii, pressure drop, and fluid viscosity of the experimental measurements using Equations |2] 
and|3](see Methods). The symbols (same as in Figure 2) correspond to those data in Figure 2 below the line / = 0, and are 
placed here by solving Equation |5] (see Methods). Surprisingly, extremely long channels (purple, red, and blue data) can 
show significant end effects. Limit of green line as / ^- oo is ^^ R/L (dashed line). 



applications of fluid-transporting nanochannels. Our results also indicate that end effect friction in bi- 
ological channel ent ryways may be a significant component in determining and controlling flow rates 



(Hummer et al 



200 ih 



Discussion: The validity of our continuum scale predictions have yet to be carefully tested for the narrowest 
channels where molecular discreteness effects may be prominent. The published flow rates from molecular 



dynam ics studies (jKalra et alll2003 



Suk and Aluru 2010 



Thomas and McGaugheyl2008 



Goldsmith and Martens 



20091) are close to the frictionless channel limit. This is in contrast to experimental reports which show some 
flow rates significantly above this limit, see Figure 2. However, the molecular dynamics studies do point 
to areas of further investigation. The p ressure drop at t he ent rance and exit of a channel are predicted by 



the continuum model to be equal, but in 



Suk and Alurul (|2010|) . are not. This may indeed be a phenomenon 



unique to molecular scale flows. However, it may also be a result of the small size of the fluid reservoirs (the 
continuum model assumes infinite reservoirs), to the unique way in which the upstream and downstream 
reservoirs are connected, to non steady-state, or even to nonhomogeneous system heating. 



Methods: The Hagen-Poiseuille equation gives the pressure drop, neglecting end effects, for fully-developed 
flow through a cylindrical channel. 



Afl 



HP 



S^LQ 



(1) 



where Q is flow rate, jU is dynamic viscosity, r is channel radius, L is channel length, and the friction factor 



/=1/^, 



where the slip length Z? quantifies th e amount by which the liquid slips along the solid channel wall 



(JMajumder et a 



2005 



Holt et al 



20061). Q diverges to infinity as / — ;■ 0. 



The pressure drop due to end effects, APg, is well approximated by the expression for flow through an 



aperture dSampsonlll 89 iHWeissbergll 19621 : ISuk and Alurull2010|) . 



AR 



s 



3Q^ 



(2) 



Using the expression for total pressure drop across a channel between two reservoirs, A/V = AP^ + APhp, 
Equations [U and [2l give a modified Hagen-Poiseuille equation, 



APt 



%^L 



3m^ 



Qt 



(3) 



where the term in brackets [. . .] is defined as the total flow resistance Rj, see Figure 3. Equation [3] can be 
rearranged to give the flow rate in the limit of large slip {b^ r). 



Q = Qq 



2 L 

i + ^t: 



(4) 



where Qq = APrr^ /3il is the flow rate through a frictionless channel (Z? — )■ oo), and the subscript T on the 
flow rate has been dropped for simplicity. Equation |4] takes into account both high slip and end effects in 
nanoscale channels, in contrast to the normal HP equation (O. By equating the flow rate without and with 
end effects (Eqs. 1 and 3), the slip length can be expressed in terms of previously-found slip lengths b- 
which neglected end effects, b = b- {1 — ^-)^^- Experimental measurements which found b- /L > 2j3n 
are unphysical. 



I 

Fabricated carbon nanotube membranes (Majumder et alll2005 



Holt et al 2006 



Hinds et al 2004 



Du et al 



201l|) are comprised of a heterogenous distribution of nanotube radii, so the total membrane flow rate is 



iV 



where A^ is the number of open nanotubes spanning the membrane and Qiji) is the flow rate 



through c 



lannel 



/■ We generate a normal distrib ution of radii, ym, based on published radii histograms (|Holt et al 



Hinds et al 



20041 : 



Du et al 



20061 : 



201l|) . From this distribution the flow rate through a membrane comprised of 



frictionless carbon nanotubes is 





liquid 


M 
(CP) 


N 


L 
(Mm) 


Gm 
(cm^/s) 


(nm) 


eM,o 

(cm^/s) 


s 


b 
(Ijm) 




water 


1 


10^ 


34 


9.62 
16.9 


<r„>=3.5 
(J.=1 


1.79 


5.37 






9.44 


- 


Majumder 
(2005) et al. 


3.410^ 


126 


12.4 


6.07 


2.04 


- 


ethanol 


1.1 


5.89 


5.06 


1.16 


- 


isopropanol 


2 


1.47 


3.10 


0.47 


24.11 




hexane 


0.3 


7.33 


20.7 


0.35 


14.66 




decane 


0.9 


0.877 


6.75 


0.13 3.99 


Majumder 
(2011) etal. 


water 


1 


2.410^ 


81 


2 10-^ 


<r„>=3.5 

o„=1.0 


8.3-10-^ 


2.40 


- 



Holt 
et al. 



2 
3 
2.8 



7.47 
2.23 
1.97 



<r„>=0.8 
o„=0.15 



4.68 



1.60 

0.48 0.58 

0.42 0.43 ^ 



Du 
et al. 



water 


1 


hexane 


0.31 


dodecane 


1.34 



0.058 
0.154 
0.082 



<r„>=5.0 
(J„=1.0 



0.112 0.52 910 _ 
0.362 0.43 630 _ 
0.084 0.97 30000 



Table 1. Data from prior work (iMaiumderet all 120051; l201ll; iHolt et all 120061; iDu et alll201 ih (left of vertical dividing 
line) and derived values (right of dividing line) used in Figure 1 . ;U is dynamic viscosity, A^ number of open nanotubes 
spanning Ic m^ of membrane, L nanotube length, and Qm total flow ra te through Icm-^ membrane. For iMaiumder et all 
(I2005L I2OIII) . Qm was computed from their published velocities For iHolt et a ll (!20 06|), Qm was computed from their 
published slip lengths. Published distributions of nanotube radii (iHolt et all 120 06; Hind s et all 120041; JDu et al|[201ll) were 
used to generate the normal distribution of radii used in the membranes, r^, with mean (r^) and standard deviation ctm, 
used in our calculations. Qm,o is the total flow rate through a 1 cm^ membrane comprised of frictionless nanotubes of radial 
distribution tm- b is the slip length from Eq. [5] Values of cS = Qui Qua > 1 do not comply with known physics of fluid 
flow. 



A' 






The slip length can then be determined from Eq. IH 



h = L 



2 ^ 
SttI-^' 



(5) 



where ^ = Qm/Qm,o is the ratio of the measured flow rate to the maximum rate that would occur through 
the same distribution of frictionless nanotubes. 



For aquaporin, data is derived from published osmotic permeabilities pf (freston et all 
in terms of hydrostatic pressure drop across the channel, 



1992|) expressed 



mpf 



where p is fluid density, kT is the thermal energy, and m is the molecular mass. The end-effect pressure 
drop is given by Eq. |2l with the channel radius chosen as the average radius through the single-file portion 
of the channel. 

Values used to produce Figures 2 and 3 and a summary of experimental measurements are given in 



Table[T] In the figures, the friction factor, /, is given in units of |^m^^ The flow resistance, Ra = IsPa/Q, 
where a = T,S for total membrane or end-effects, respectively, and APa is the pressure drop. 

In Figure 1, the rate of energy dissipation was calculated from a numerical simulation of the Navier- 
Stokes equations using an axially symmetric geometry of a cylindrical channel between two large reservoirs 
with a flow rate of 10-^|J.m-^/s. A perceptually-balanced color scale in MATLAB was chosen for visual accu- 
racy to reduce luminance variation artifacts of the standard rainbow colormap. In the simulations, changing 
from slip to no- slip boundary conditions on the reservoir walls had negligible effects on flow rates. Addi- 
tionally, interactions between neighboring channel openings were investigated by adding a second nearby 



channel. At the separation for adjacent doub 



such as may occur in fabricated membr anes (IHolt et al 



agreement with theoretical predictions dWang 



e-walled carbon nanotubes and internal channel radius 1 nm. 



1994) . 



20061) . no significant interactions were observed, in 
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